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^ ' Abstract 

o 



Very recently a new algorithm to the nonnegative single-source shortest 
lO ■ path problem on road networks has been discovered. It is very cache- 

efficient, but only on static road networks. We show how to augment it 
to the time-dependent scenario. The advantage if the new approach is 
that it settles nodes, even for a profile query, by scanning all downward 
edges. We improve the scanning of the downward edges with techniques 
C/2 ' developed for time-dependent many-to-many computations. 

o ' 



1 Introduction 

> 

Q\ ' The new algorithm for the nonnegative single-source shortest path problem on 

■ road networks [5] uses contraction hierarchies [8]. First, it performs a forward 

upward search, and then it processes all nodes in descending order of importance. 
The algorithm is so efficient because it can order the edges very cache-efficient 
and can be parallelized. We consider the time-dependent scenario, where the 
travel time depends on the departure time, the edge weights are complex travel 
time functions (TTFs). Such a TTF maps a departure time to the travel time. 
We want to solve the problem of computing the travel time profiles from one 
source node to all other nodes in the graph. A travel time profile is a TTF 
, that maps the departure time at the source to the earliest arrival time at the 

5-h ' target node. Significantly the most work on solving our problem is to process 

the TTFs, the traversal of the graph takes negligible time. So the cache-efficient 
order of the edges brings is no longer a significant speed-up. But the way the 
edges are processed allows to prune a lot of expensive TTF operations by using 
approximate TTFs. 

1.1 Related Work 

Many new algorithms for the time-dependent point-to-point shortest path prob- 
lem have been developed recently. We refer to [6] for an overview. Also, a 
time-dependent version of contraction hierarchies (TCH) [1, 2] exists, and we 
use it for our new algorithm. It was augmented to compute travel time tables 

by [7]. 
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2 Preliminaries 



2.1 Time-Dependent Road Networks 

Let G = (V,E) be a directed graph representing a road network. 1 Each edge 
(u, v) G E has a function / : R — > R>o assigned as edge weight. This function / 
specifies the time /(r) needed to reach v from u via edge [u, v) when starting at 
departure time r. So the edge weights are called travel time functions (TTFs). 

In road networks we usually do not arrive earlier when we start later. So 
all TTFs / fulfill the FIFO-property: W > r : t' + f(r') > r + /(r). In this 
work all TTFs are sequences of points representing piecewise linear functions. 2 
With |/| we denote the complexity (i.e., the number of points) of /. We define 
/ ~ g Vr : /(r) - 5 (r) for ~G {<, >, <, >}. 

For TTFs we need the following three operations: 

• Evaluation. Given a TTF / and a departure time r we want to compute 
/(t). Using a bucket structure this runs in constant average time. 

• Linking. Given two adjacent edges (u,v), (v,w) with TTFs /, g we want 
to compute the TTF of the whole path (u ->/ v — > g w) . This is the TTF 
g*f : t n- <7(/(t) + t) + /(t) (meaning g "after" /). It can be computed in 
0(|/| + \g\) time and \g* f\ G 0(|/| + \g\) holds. Linking is an associative 
operation, i.e., / * (g * h) = (/ * g) * h for TTFs /, g, h. 

• Minimum. Given two parallel edges e, e' from u to v with TTFs /, /', 
we want to merge these edges into one while preserving all shortest paths. 
The resulting single edge e" from u to v gets the TTF min(/, /') defined 
by r H min{/(r), /'(r)}. It can be computed in 0(|/| + |/'|) time and 
|min(/,/')| G 0(|/| + |/'|) holds. 

In a time-dependent road network, shortest paths depend on the departure 
time. For given start node s and destination node t there might be different 
shortest paths for different departure times. The minimal travel times from s 
to t for all departure times r are called the travel time profile from s to t and 
are represented by a TTF. 

Approximations. Give a TTF /. A lower bound is a TTF with /^ < / 
and a lower £ -bound if further (1 — e)f < f^. An upper bound is a TTF f' 
with / < f' and an upper e -bound if further f> < (1 + e)/. An e- approximation 
is a TTF ft with (1 - e)f < f < (1 + e)f. Approximate TTFs usually 
have fewer points and are therefore faster to process and require less memory. 
To compute £-bounds and ^approximations from an exact TTF / we use the 
efficient geometric algorithm described by Imai and Iri [9] . It yields a TTF with 
minimal number of points for e in time 0(|/|). 

2.2 Time-Dependent Profile Dijkstra 

To compute the travel time profile from a source node s to all other nodes, 
we can use a label correcting modification of Dijkstra's algorithm [10]. The 
modifications are as follows: 

1 Nodcs represent junctions and edges represent road segments. 

2 Here, all TTFs have period il = 24h. However, using non-periodic TTFs makes no real 
difference. Of course, covering more than 24h will increase the memory usage. 
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• Node labels. Each node v has a tentative TTF from s to v. 

• Priority queue (PQ). The keys used are the global minima of the labels. 
Reinserts into the PQ are possible and happen {label correcting). 

• Edge Relaxation. Consider the relaxation of an edge (u, v) with TTF /„„. 
Let the label of node u be the TTF /„. The label /„ of the node v is 
updated by computing the minimum TTF of f v and f uv * f u . 

2.3 Time-Dependent Contraction Hierarchies 

Hierarchies. In a time- dependent contraction hierarchy [1] all nodes of G are 
ordered by increasing 'importance' [8]. In order to simplify matters, we identify 
each node with its importance level, i.e. V = l..n. 

Now, the TCH is constructed by contracting the nodes in the above order. 
Contracting a node v means removing v from the graph without changing short- 
est path distances between the remaining (more important) nodes. This way we 
construct the next higher level of the hierarchy from the current one. A trivial 
way to contract a node v is to introduce a shortcut edge (w, w) with TTF g * / 
for every path u — v — > g w with v < u,w. But in order to keep the graph 
sparse, we can try to avoid a shortcut (tt, w) by finding a witness - a travel 
time profile W from u to v fulfilling W < g * / . Such a witness proves that the 
shortcut is never needed. The node ordering and the construction of the TCH 
are performed offline in a precomputation and are only required once per graph 
independent of S and T. 

Queries. In the point-to-point scenario, we compute the travel time profile be- 
tween source s and target t by performing a bidirectional time-dependent profile 
search in the TCH. The special restriction on a TCH search is that it only goes 
upward, i.e. we only relax edges where the other node is more important. This 
property is reflected in the upward graph G-|-:= (V, E^) with E^:= {(it, v) € E | u < v} 
and, the downward graph G^:= (V,E±) with E±:= {(u, v) G E \ u > v}. Both 
search scopes meet at candidate nodes u giving lower/upper bounds on the travel 
time between source and target, allowing us to prune the following profile search. 
The bidirectional profile search computes forward TTF /„ and backward TTF 
g u representing a TTF g u * f u from source to target (though not necessarily an 
optimal one). The travel time profile is mhi{g u * f u \ u candidate}. 

3 Our Algorithm 

Given a source node s, we want to compute for each node u in the graph the 
travel time profile S(u) from s to u. We initialize S(s) := (r M> 0) and all other 
TTFs 6(u) = (t4 oo). Then, we perform a forward search from s in Gf and 
update the tentative travel time profile S(u) for all nodes visited by this search. 
Now, we process all nodes in the graph by descending importance level and 
compute 

S(u) :— min (S(u), min {f v * S(v) \ v — >f v «£ E±} ) (1) 

So essentially, we scan through all incoming downward edges v — >f u of u, 
link them to S(v) and build the minimum. Doing this naively is very costly, as 
the travel time functions arc complex. As not all of the nodes v contribute to the 
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minimum in the end, we improve the performance by using pruning techniques 
developed for the time-dependent travel time table- computation. Assume that 
we have lower/upper e-bounds /j / /J for /„ and S(v)^ / S(vy for S(v). Then 
we can use Algorithm 1 to accelerate the computation of 6(u). 

Algorithm 1: BuildMinimum(u) 

1 5 := min v ^. U £E, {max/„ + maxi5(w)}; // upper bound based on 
maxima 

2 {v. ~ ► ') := ar g mm u^ f ueE l { mm / + min<5(u)}; // minimum node 

3 5^ := fy * 5(v)^; II upper bound based on approximate TTFs 

4 5 := mm(6, max^); // tighten upper bound 

5 foreach v —>f v U G E± do // loop over all downward edges 
if min/t, + mia6(v)^ < 5 then // prune using minima 

0^ := min (<r, /J * ; // update upper bound 

8 S(u) := min (S(u), f v * 6(v)) ; II tentative travel time profile 

9 foreach v —tf v u G E± do // loop over all downward edges 
if * 8{up > S') then // prune using lower bounds 

8(u) := min(5(u), f v * 8(v)); II update travel time profile 

We pass three times through the incoming downward edges v — Vf v u G E±. 

1. In Line 1 we compute an upper bound 6 based on the maxima of f v and 
S(v). Also, in Line 2 we compute the edge with minimum sum of the 
minima of /„ and S(v). This edge is usually very important and a good 
starting point to obtain an tight lower bound. 

2. In Lines 3-7 we compute an upper bound 5^ based on the upper e-bounds. 
This bound is tighter than the one based on the maxima. 

3. In Lines 8-11 we compute the travel time profile and use the upper bound 
5' for pruning. So we only execute the very expensive link and minimum 
operations on /„ and 5(v) at Line 11. 

In comparison, a Profile Dijkstra would update S(u) gradually when he pro- 
cesses node v, and thus cannot pass through all the downward edges several 
times to compute the intermediate upper bounds S and 5*. 
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4 Improvements 

Some of the improvements from [5] can be applied, especially the reordering of 
nodes, and the parallel processing. However, we did not adopt SIMD instruc- 
tions or GPU, as we now operate on complex TTFs. 



5 Core-based Computation 

An interesting observation is that [5] cannot be used like Dijkstra's algorithm to 
compute the distances to close-by nodes, as the nodes are no longer processed 
in order of increasing distance. However, we can compute the distances only for 
a core of the contraction hierarchy, that is a number of k most important nodes. 
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algorithm 



threads prune e [%] 



query time [s] 



Dijkstra 



116 



TCH 
TCH 
TCH 
TCH 
TCH 



1 
1 
1 
1 
1 



0.1 
0.01 



10 

1 



105 
108 
48.2 
32.6 
35.6 



TCH 
TCH 
TCH 
TCH 



2 
4 
6 

8 



0.1 
0.1 
0.1 
0.1 



16.0 

8.8 
6.5 
5.7 



Table 1: Performance 



Computing only the distances to all core nodes is faster, especially in the 
time-dependent scenario. Also, it takes much less space, as the TTFs usually 
contain thousands of points. 

An application of the core-based computation is the computation of arc flags 
[5, §7.2], but now only for a core. In [3] it is observed that this brings significant 
speed-ups for the static scenario. However, in the time-dependent scenario, 
upper/lower bounds used to compute exact arc flags, as using exact TTFs is 
too time-consuming [4]. These arc flags are very weak, and simple heuristic 
computation of the arc flags using time-sampling significantly accelerate the 
query, but provide no approximation guarantee [4] . So we expect that when we 
are able to compute exact arc flags with exact TTFs, this provides both a fast 
and exact query. 

6 Experiments 

Input. We use a real- world time-dependent road network of Germany with 4.7 
million nodes and 10.8 million edges, provided by PTV AG for scientific use. 
It reflects the midweek (Tuesday till Thursday) traffic collected from historical 
data, i.e., a high traffic scenario with about 8 % time dependent edges. 

Hardware/Software. The experiments were done on a machine 3 with two 
Intel Xeon E5345 processors (Quad-Core) clocked at 2.33 GHz with 16 GiB of 
RAM and 2x8 MiB of Cache running SUSE Linux 11.1. We used the GCC 4.3.2 
compiler with optimization level 3. 

Basic setup. We use a preprocessed TCH as input file [1]. We use a core-size 
of 10 000, and select 100 core nodes uniformly at random as source of our query. 

The experimental results are presented in Table 1. The number of threads 
corresponds to the number used to answer a single query. 

3 The machine used in [2, 7] is currently repaired, we plan to publish results for this machine 
later. 
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7 Conclusion 



Wc presented an efficient algorithm for time-dependent one-to-all computation. 
By applying it to a core, it can be used to accelerate precomputation of speed-up 
techniques. 
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